Constraints and tests of the OPERA superluminal neutrinos 
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The superluminal neutrinos detected by OPERA indicates Lorentz invariance violation (LIV) of 
the neutrino sector at the order of 10~^. We study the implications of the result in this work. 
We find that such a large LIV implied by OPERA data will make the neutrino production process 
TT — i- /.t + z^^ kinematically forbidden for neutrino energy greater than about 5 GeV. The OPERA 
detection of neutrinos at 40 GeV can constrain the LIV parameter to be smaller than 3 x 10~^. 
Furthermore the neutrino decay in the LIV framework will modify the neutrino spectrum greatly. 
The atmospheric neutrino spectrum measured by IceCube can constrain the LIV parameter to the 
level of 10~^^. The future detection of astrophysical neutrinos of Galactic sources is expected to be 
able to give even stronger constraint on the LIV parameter of neutrinos. 



Introduction — Recently the OPERA collaboration re- 
ported an astonishing discovery that the neutrinos from 
CERN to the Gran Sasso Laboratory may be superlu- 
minal The arrival time of the muon neutrinos with 
average energy ^17 GeV is earlier by 



St = (60.7 ± 6.9(stat.) ± 7.4(sys.)) ns 



(1) 



than that of light, which corresponds to a difference of 
the muon neutrino velocity and the speed of light of 

6^^ = {v^~c)/c = (2.48±0.28(stat.)±0.30(sys.))xl0~^ , 

(2) 

with significance level of 6a. The energy dependence 
study of OPERA shows that the velocity difference Siy-y 
is almost independent of energy. It was found that, by 
splitting the events into two groups with energies be- 
low or above 20 GeV, the velocity difference is given by 
6^^ = (2.16 ± 0.76 ± 0.30) x lO'^ for (E) = 13.9 GeV 
and 6^^ = (2.74 ± 0.74 ± 0.30) x 10"^ for {E) = 42.9 
GeV. The results above have been confirmed recently by 
a test performed using a beam with a short-bunch time- 
structure allowing to measure the neutrino time of fiight 
at the single interaction level. The new analysis show 
consistent result with that at the first version. 

Actually the earlier experiments also gave similar re- 
sults about the velocity of neutrinos, but with lower sig- 
nificance. In 2007 the MINOS experiment at FERMT 
LAB measured the neutrino speed for energies around 3 
GeV and found 



<5^^ = (5.1 ±2.9) X 10" 



(3) 



which seems to be consistent with that by OPERA. 
Results from even earlier experiments at FERMILAB 
showed a similar value of 6^^ at the order of 10~^ 3. 
However, the supernova neutrinos from SN1987A set a 
much more stringent limit of < 2 x 10~^ for tens of 
MeV electron neutrinos [J]. 

Soon after the OPERA result was reported there were 
many works to discuss the implications and possible ex- 
planations to it 041] ■ The analysis in these works seem 
to get a consensus that the effect should be neutrino fla- 
vor independent so that the coherent neutrino oscillation 



is not destroyed. Furthermore, the effect as a function 
of energy should be a step-like function in order that it 
is large and energy independent at OPERA and FER- 
MILAB energies (>GeV), while at low energies it is very 
small to be consistent with the SN1987A result. 

Phenomenologically such a superluminal neutrino can 
be easily accommodated in a Lorentz invariance viola- 
tion (LIV) model. Coleman and Glashow (hereafter we 
refer as CG) constructed such a model where LIV was 
expressed as a modified energy-momentum relation 



(4) 



assuming that the Lorentz invariance is violated pertur- 
bati yely in the context of conventional quantum field the- 
ory [10| . In such framework, the neutrino velocity is given 
by Wj, = 1 — m'^/2\p\'^ +^/2 taking c = 1. For high energy 
neutrinos we have (5i/-y = ^/2. 

In this work we will assume that the LIV term is fla- 
vor independent and given by the form in Eq. Q at high 
energies. At low energy, the LIV term is suppressed by 
some specific models, such as that given in [3, Q, to be 
consistent with SN1987A data. As studied in the CG 
paper, many phenomenologies changed significantly. Es- 
pecially for a large ^ of the order 0(10"^), the processes 
7r+ — > /i+ + i/fj^ and fi^ — > + 1^^ + />e are even forbidden 
so that no high energy neutrinos could be produced. In 
the next section wc will show the constraints on the LIV 
parameters by such processes. 

In addition, Eq. (jj]) can also induce the neutrino ra- 
diative decay — > I'j + "f and the three body decay 
i^i — > I'i + + e or i/i Ui + e + e. Such processes 
will lead to the fast energy loss of neutrinos once the 
threshold is reached. We will discuss the consequence on 
the atmospheric and astrophysical neutrino spectra when 
considering such processes. 

Constraints from neutrino related processes — We first 
consider the process 7r+ /i+ + ly^. For ^ > 0, the LIV 
term in Eq. ^ is equivalent to an "effective" mass term. 
If this effective mass is too large, we will expect kinemat- 
ics forbiddance of the pion decay. To derive the threshold 
condition wc get the effective mass of superluminal ncu- 



2 



trino following CG [T^l as 



(5) 



where we take the LIV parameters of pion and muon 
to be zero as they can be much smaller than the neu- 
trino LIV parameter. We will give a discussion on the 
LIV parameters of pion and muon in the next section. 
Neglecting neutrino mass rrii, we get the threshold con- 



dition for 7r+ decay 7r+ fi'^ 



off 



y^7(l + C) • the laboratory frame we have 



< (m^-m^)v/r + T7l 



(6) 



This inequality means that neutrinos generated by pion 
decay have a maximum energy beyond which the decay 



process tt^ 



can not occur at all. 



Similarly the processes K'^ — ^ fj,~^ + 1'^ , M ~^ f + e -I- 
and T ^ Vr + e(/^) -I- j^e(/i) should also be forbidden at 
high energies. In Fig. [T]we show the constraints on the 
neutrino velocity {vi, — c)/c for different neutrino energies 
with different processes. The region above the line should 
be excluded if we detect the neutrinos with correspond- 
ing energies. The experimental data by FERMILAB79, 
MINOS and OPERA are also shown. From Fig. [1] we 
can see that the process tt — >■ fii'p,, through which the 
OPERA neutrinos are produced, gives the most stringent 



constraint. Process fi ^ i'^ 



gives a similar but 



a little weaker constraint. It is shown that the detection 
of neutrinos with energies as high as 43 GeV by OPERA 
can set a constraint < 3 x 10~^, while the ~ 200 GeV 
neutrinos detected at Fermilab give 5^^ < 2 x 10~*. 
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FIG. 1: The neutrino velocity constrained by the neutrino 
energy, for processes -k{K)^ — > fJ-iT) — >■ i'^(t) + e(M) + 



Actually the highest energy neutrinos detected are the 
atmospheric neutrinos at IceCube, which reach ^ 400 
TeV These neutrinos are generated by the same 

process as OPERA neutrinos, i.e. through — >■ tt{K) — >■ 
jJL+v^. From Fig. [l]we can see that the IceCube neutrinos 
can constrain LIV parameter down to ~ 4 x 10~^'^. 

In the following we extend the CG model with a more 
general case that the dispersion relation is modified as 



E' 



m2 + |pl2+/3|pl" 



(7) 



We use the same particle decay processes to constrain the 
model parameters. The constraints from neutrino energy 



for 7r+ 



are shown in Fig. [21 for different values 



of a. The results show that larger a will give weaker 
constraints on the parameters. 



10-^ 
10"^ 
10"° 
10"^ 
10"' 
10"^ 
10"' 



a=6 
ra=10 




IceCube | 
atm neutrino 



10" 10' 10' 
Ev (GeV) 



10' 1 0'' 10° 



FIG. 2: The neutrino velocity constrained by the neutrino 
energy for the process tt^ — >■ /i^z^^, for different values of a. 



The atmospheric and astrophysical neutrino spectra — 
We have discussed the constraints on the LIV parameter 
of neutrinos according to the processes of tt — > /i + 
and ^ + e + Ve- These processes can only constrain 
the difference of the LIV parameters between the parent 
particle and the daughter particles. Therefore to give 
stronger constraint on the neutrino LIV parameter we 
have to give constraints on tt and /-t first. 

For superluminal 7r+ it will lose energy quickly by 
radiative emission tt^ — > tt'^j. IceCube experiment 
has measured the atmospheric neutrino spectrum up to 
~ 400 TeV, which fits very well with the model calcu- 
lation [lT|. This means that the high energy charged 
pions do not lose much energy before they decay to 
neutrinos. This can give a strong constraint on the 
process 7r+ — > 7r+7. The threshold of this process is 
i?7r > TO^/ -^/cj — = m-^/y/^ [l^l . This leads to an 
upper bound of < 10^^** when the maximum energy 
of pion is about 2 PeV. 

If the maximum attainable velocity of muon is smaller 
than c, it is still possible to have the process tt — > -|- 
with superluminal Vf_i. In such a case the LIV parameter 
is constrained by the process 7 — >■ As the 7 e nerg y as 
high as 80 TeV has been detected at Crab Nebula [i2| , "we 
have l^pl < I0~^^. Therefore the LIV parameter of neu- 
trinos can be constrained at most to the order ~ 10~^^ 
according to pion decay [isj . More stringent constraint 
on the flavor dependent LIV parameter was set in [l^ . 

In the following we study how to give more stringent 
constraints on the neutrino LIV parameter. For superlu- 
minal neutrino with the LIV dispersion relation Eq. (U), 



it will lose energy through radiative decay Vi 



+7 



or the three body decay. The three body decay can be a 
charged current (CC) process i^,, — > /i + + e or a neu- 
tral current (NC) process + e + e. Among these 



3 



three processes, the CC process has the largest decay 
rate, but the threshold energy is also the highest. The 
radiative process has the lowest decay rate with lowest 
threshold energy. In the OPERA neutrino energy loss 
is studied when propagating from CERN to Gran Sasso 
taking ^ = 5 x 10^^. As we have shown that neutrinos 
with energy up to tens of GeV can not be produced if ^ 
is so large, we will study the energy loss process with a 
much smaller ^. As argued above if muons are involved, 
the LIV parameter of neutrinos can only be constrained 
to the order 10~^^. Therefore in the following we focus 
on the NC process which involves electrons in the final 
states, as discussed in @. The processes e — >■ 67 and 
7 — > ee give |^e| < 10""'^^ by observation of the Crab 
Nebula fia. 
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FIG. 3: The expected atmospheric neutrino spectra with LIV 
parameters ^ = 10"^°, 10"" and 8 x 10"" respectively. The 
measured data by IceCube is also plotted [lH . 



The calculation of the three body decay rate is very 
complicated if the LIV effect of final state particles is 
taken into account. Working out the integration of final 
state phase space, we get 



energy interval dilation dEo/dE is 

r*=^/^ db[Eo{t')] 



dEo 
dE 



exp 



dEo{t') 



dt' 



(10) 



where b{E) = —dE/dt is the energy loss rate in Eq. ([8|). 
Note here we consider the local neutrino sources. For 
cosmological sources as discussed in [l^ Eqs. (9) and 
(10) are more complicated. 

We first study the atmospheric neutrinos. In [ll| the 
up-going neutrinos from 100 GeV to 400 TeV were de- 
tected by IceCube. We take a rough travel distance of 
lO'* km for the neutrinos with zenith angle 1 24° — 180°. In 
Fig. [3]we show the neutrino spectra with LIV parameters 
^ = lO"^'', 10^^^, and 8 x 10^^"^ respectively. It is shown 
that for specified distance of the source, there is a maxi- 
mum energy of the neutrinos above which the energy loss 
is too fast that the neutrinos can not reach us with high 
energies. According to the energy loss rate in Eq. ([8]) , we 
have EgcV = (^o:Lv + 5 x lO^^s^scc)"'/^ For Eo ^ 00, 
£^max oc t"!/^ (d/c)"i/^ For ^ = 8 X 10"" we have 
Emax ~ 400 TcV for atmospheric neutrinos. According 
to the fact that IceCube has detected atmospheric neu- 
trinos up to 400 TeV, we can set a bound of the LIV 
parameter of neutrino ^ < 8 x 10" 
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(8) 

where (Ai?) and F arc the averaged energy loss of each 
decay and the decay width. A detailed calculation in the 
laboratory system has been given in [isl ] . 

To calculate the neutrino spectra with the energy loss 
due to decay, we follow the way of calculating the spec- 
tra of ultra high energy (UHE) protons propagating in 
the cosmic microwave background (CMB) field [i6|. The 
detected neutrino spectra can be written as 



<i>{E) 



qoiEg) dEp 

47rd2 dE ' 



(9) 



where 50(^0) is the production spectra per unit time at 
the source, d is the distance of the source, Eq and E are 
the energies at the source and detector respectively. The 



FIG. 4: Expected neutrino spectra with LIV parameter ^ = 
10" ^^ for sources at distances 0.3, 1.0 and 3.0 kpc respec- 
tively. The source spectrum is assumed as E~^. 



Then we study the effect of energy loss on the as- 
trophysical neutrino spectrum. Because these neutrinos 
will propagate very long distance before reaching the 
Earth we expect the spectral measurement can probe 
even smaller LIV parameters of neutrinos. Up to now, no 
high energy astrophysical neutrinos were detected. How- 
ever, as the well believed sources of Galactic cosmic rays 
(CRs), supernova remnants (SNRs) are expected to be 
good candidates of astrophysical neutrino sources. Fig. 
m gives an illustration of the astrophysical neutrino spec- 
tra taking into account the LIV effect with much smaller 
parameter ^ — 10"^^. Three distances of the sources, 
0.3, 1.0 and 3.0 kpc are adopted. In the calculation the 
source spectrum is assumed to be cx E~^. We can see 
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that for sources with distance '-^kpc the cutoff energy is 
about several tens to 100 TeV for this small LIV param- 
eter. The 7-ray observations of the Galactic SNRs show 
that the accelerated energies of particles can well exceed 
TeV, which enables the production of neutrinos with en- 
ergies much higher than TeV. Therefore the observation 
of cutoff at neutrino spectra in the future can set a very 
strong constraints on the LIV parameter. 

Finally the Greisen-Zatsepin-Kuzmin (GZK) neutri- 
nos are expected when the ultra-high energy cosmic rays 
(UHECRs) interact with the CMB photons and generate 



the GZK cutoff of UHECRs [17[. It is expected that the 



IceCube might discover such neutrinos in the near fu- 
ture. However, as shown above the decay mean free path 
for tens TeV neutrinos is only of kpc scale. We might 
have no chance to detect any neutrinos from extragalac- 
tic sources if they are superluminal. LIV at the proton 
sector may also suppress the GZK neutrinos jl8| . 

Conclusion — We studied the implications of the claim 
of superluminal neutrinos by the OPERA collaboration. 
The observation generally indicates LIV in the neutrino 
sector. In such a framework, we find the neutrino pro- 
duction process 7r+ — >■ fJ.'^i'^ becomes kinematically for- 
bidden for > 5 GeV if taking S^^ = 2.5 x 10~^ as 
shown by OPERA data. The detection of ^ 40 GcV 
and ^ 200 GeV neutrinos by OPERA and Fcrmilab can 
constrain the LIV parameters to 3 x 10~^ and 2 x 10~® 
respectively. 



Then we discuss if it is possible to constrain the neu- 
trino LIV parameters by studying the three body decay 
process of neutrinos. Such process will modify the at- 
mospheric and astrophysical neutrino spectra. Because 
there is no spectral cutoff of the atmospheric neutrinos up 
to ^ 400 TeV as measured by IceCube [lH, the neutrino 
LIV parameter can be constrained to less than 10""'^^. 

We further studied the astrophysical neutrino spectra 
from distant sources, for example the Galactic SNRs, 
which arc usually taken as the sources of Galactic cos- 
mic rays. The long travel distance makes it possible to 
probe even smaller LIV parameters. As an illustration, 
we find that for ^ ~ 10~^^ and source distance ~ kpc, 
the decay process is expected to result in a distinctive 
spectral cutoff of the neutrino spectra at energies about 
several tens TcV. This effect can be tested through the 
future detection of the neutrino sources. 
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